Effective cross-Kerr nonlinearity and robust phase gates with trapped ions 
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We derive an effective Hamiltonian that describes a cross-Kerr type interaction in a system involving a two- 
level trapped ion coupled to the quantized field inside a cavity. We assume a large detuning between the ion and 
field (dispersive limit) and this results in an interaction Hamiltonian involving the product of the (bosonic) ionic 
vibrational motion and field number operators. We also demonstrate the feasibility of operation of a phase gate 
based on our hamiltonian. The gate is insensitive to spontaneous emission, an important feature for the practical 
implementation of quantum computing. 
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Nonlinear effects are specially relevant to quantum optics as well as to applications in quantum information and computation. 
In particular, the cross-Kerr coupling is at the basis of various schemes, such as quantum non-demoUtion measurements {1101; a 
Fock state synthesizer jsl, entanglement purification |4||, teleportation 0], quantum codewords 0], a proposal of realization of 
a CNOT gate Q], and Bell-state detection 0|, for instance. Large cross-Kerr nonlinearities are usually required for the actual 
implementation of those schemes, but the values of the Kerr coefficients in the available nonlinear media are in general relatively 
small, although some progress based on atomic gases has been achieved . It is therefore of interest to seek alternative 

ways of obtaining nonlinear interactions of that type. The cross-Kerr interaction may involve the effective coupling of two modes 
of the field Hill , or even two sub-systems of different nature, such as an atom and a field Hill . 

Thinking about specific applications, and given the extraordinary development of the field of quantum computing, it is im- 
portant to identify physical systems and situations that could be suitable for the realization of efficient quantum logic gates. A 
rather interesting physical system comprises a single trapped ion coupled to the quantized field iHsll . That configuration forms 
a tripartite quantum system: the internal electronic degrees of freedom; the quantized center-of-mass motion; and the quantum 
cavity field. Several aspects of that system connected to quantum information and computation have been already investigated, 
e.g., a scheme for generation of entangled states involving the ionic vibrational motion and the quantized field lll4ll . as well as 
an alternative proposition of implementation of a CNOT gate fisll . As a matter of fact, the possibility of coherent manipula- 
tion of trapped ions (qubits suitable for storing quantum information) coupled to photons (qubits suitable for carrying quantum 
information) iHfill opens new possibilities for quantum information processing fl7|. 

In this paper we derive, from a general hamiltonian involving the oscillating ion interacting with the quantized field, an 
effective Hamiltonian which couples the quantized field and the ionic motion in a way similar to what it is found in nonlinear 
optics. More specifically, we obtain a coupling which is basically the product of bosonic number operators referring to the 
vibrational motion (0^0) and the quantized field (&^fe), yielding an interaction hamiltonian of the form Hint ~ Aa^dfotfo. In 
the second part of our paper, we suggest an application of our Hamiltonian for quantum computing, presenting a proposal of 
implementation of a quantum phase gate. In our scheme, the ion is initially prepared in its ground (electronic) state \g), remaining 
in that state during the gate operation, which means that the gate will be quite robust against atomic spontaneous emission. 

We consider a two level ion (with atomic frequency cja) trapped in a harmonic potential (frequency v) inside a high finesse 
cavity, and coupled to a single mode (frequency uJc) of the quantized cavity field. The ion has its motion well cooled down and 
therefore it should described quantum mechanically. The Hamiltonian for that system reads 11311 

H = ua^ a + ojcb'^b + — WaO'z + g{<^+b + b^o'-) cos7^(d^ + d), (1) 

where d^^ (d) are the creation(annihilation) operators relative to the excitations of the center-of-mass oscillatory motion, fet (5) are 
the creation (annihilation) operators of photons in the field mode, and g is the ion-field coupling constant. The operators a^, cr+ 
((T_) are the population difference and excitation (de-excitation) atomic operators, respectively. In order to derive a Kerr-like 
interaction we need to perform a couple of approximations. By one side, we note that the factor cos ?y(dt + d) in Hamiltonian 
([0 may be expanded in powers of the operators d and d^ giving then rise to different processes changing either energy or phase 
of the center-of-mass motion. It is then mandatory to get rid of the energy changing terms, keeping just terms proportional to 
powers of the number operator d^d. In the interaction picture, the Hamiltonian Q is given by 

Hi = g{a+be'^* + a A^e-'^')cosT]{a'fe"''^ (2) 
where we defined A = cjg, — ujc- Making use of the relation e^+^ = e~^^'^^^'^e^e^ , valid for [A, B] = 7, where 7 is an 
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arbitrary complex number, we have 

Hi = ga+bY,F{a\a- a, /?)e*[A+-("-/3)]* + ^^_5t ^ j^(^t^ a; a, /3)e-'[^+'^("-«]*, (3) 

a,!i a,l3 

where we defined the function 

F{a\a;a,P) = + {-^r,r+f']a^'^a^ . (4) 

Now, we perform the first approximation. By analyzing the temporal dependence in Hamiltonian (|3}, we may carefully choose 
the frequencies of the system in order to perform a suitable rotating wave approximation. In this approximation, rapidly oscil- 
lating terms are dropped out and just the slow ones are kept. The terms in oscillate in time either as e^*'^* or Q^M^+'''')t^ 
with k integer In the regime A <C and A ^ kv, the terms oscillating with lower frequencies are those proportional to e^"^*. 
This occurs if a = /3, or 

F{a\a;a,a)^e '—^ , (5) 

ar 

and then 

H = va)a + ujcb^b + + 5(^+5 + a ^b^ ) f [a'' a) , (6) 

with 

f(a^a) = e-'''/^ . j^^^) ,^ (7) 

where : Join) : is just the normally ordered zeroth order Bessel function of the first kind 1 18], with n = 2r]VdJ&. We note 
the presence of b and 6^ which must be replaced by terms proportional to wb. This may be accomplished by considering the 
large detuning limit A 3> g. We proceed and write down the Heisenberg equations of motion for the atomic and field transition 
operators. From ([Q one obtains 

*~37" = -i^acr+ + gb'^^zfia''a) (8) 
at 

= ujj) + ga^fia)a) (9) 
= ,ya + gia+b + a^b^)^fia^a). (10) 
It is now suitable to move to a new frame using the transformation 

(11) 
(12) 
(13) 

so that the above equations of motion can be rewritten as 

= gPaJia^a)e-^^' (14) 

^^ = ga^fia^a)e-^^' (15) 
at 

= .g(^+6e'^* + a„5te-^*)e— *A/(ata). (16) 
at oa' 

In the large detuning limit, i.e. for A = cj^ — cjf, 3> g, we may perform the second order approximation which consists in 
assuming that the operators vary slowly when compared to the relevant frequencies involved. This allows us to integrate il4l 
obtaining 

a+ = ^-^b^aJia^a). (17) 
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After substituting J17> into ( I15t . it results 



4 ^ ^-rfi^^^Yba.- (18) 




This equation of motion can be obtained from the effective Hamiltonian 

2 

H^s = i^a^a + uj^pb + ^(7^ + ^f{a)a){a+a^W - a^a+b^'b). (19) 

One could wonder now what is the validity of such effective Hamiltonian. Its limit of application may be found by using 
ordinary time-dependent perturbation theory. Once the effective Hamiltonian does not allow any transition between bare states 
we want the transitions \g, ric, rriv) |e, (n — l)c, m^) to be unlikely to happen. Mathematically, the transition probability is 
given by 

Pit) = |(e, (n - l)c, ml\Uit)\g, n„ m^) (20) 
where U (t) is the time evolution operator In first order, 

U{t) = 1-1 [ H{t)dt, (21) 
Jo 

being H{t) the original Hamiltonian ([l) in the interaction picture. After a lengthy calculation, the wanted transition probability 
is found to be 

(22) 

where A(t) and B(t) are oscillating functions of time. It is then clear that the effective Hamiltonian is to be safely applied when 
the conditions A kjy (fc = 0, ±1, ±2, . . .) and g <C A are fulfilled, as well as A ^ (for m). 

The Hamiltonian in eq. M9\ has novel and interesting features. It involves products of the number field operator h'^h with the 
vibrational number operator a^a, within the Bessel function Jo, characterizing a cross-Kerr type interaction between two bosons. 
The Stark shifts resulting from the large detuning condition we have assumed, gave rise to the interaction terms characteristic 
of the cross-Kerr interaction. The key point here is that the vibrational number operator a^a comes out multiplying the light 
shifts. For a sufficiently small value of the Lamb-Dicke parameter -q (Lamb-Dicke regime), we may truncate the Bessel function 
in order to obtain a term proportional to rj^a^a. From Q we note that 

f{a)a)^l-'q^-2ifa^a, (23) 

and then, if the ion is initially prepared (internal levels) in its ground state we may rewrite M9\ in that approximation as 

H^^ = va)a + CjJ)^b + \a)aPb, (24) 

where Cjc = + ??^3^/A — g^/ A is a shift in the field frequency and the coupling constant between the bosonic modes is 
A = 2rf'g^ / A. The interaction term of ( I24> is precisely the sought Kerr-type interaction in its well known form. By careful 
choice of the 77, g and A it is possible to obtain a considerable non-linearity. According to fl^ . a realistic choice for the ion-field 
coupling could be g = 2tt x 1.51MHz. In the Lamb-Dicke regime, we may consider 77 = 0.05 and for the high detuning limit 
A = IO5' is also adequate. These values lead to A « 5 kHz. 

Now, we would like to show the phase gate operation using our effective hamiltonian. We would like to point out that other 
few proposals may be found in the literature, involving trapped ions cavity QED lEHl or optical lattices |22]. There are 

also some experimental realizations 1E3I1 . Here the quantum subsystems involved are the ionic center-of-mass (phonons) and the 
quantized field (photons). Considering that the ion is in its electronic groud state, the Hamiltonian J19> assumes the following 
form in the interaction representation 

H,s = -^f{a)a)Kb, (25) 

whichis adequate for producing general phase shifts. Because [a^OjiJeff] = [i>^b,Hcfi\ = 0, we may write the evolution operator 
U{t) = exp{-iHcst/fi) in diagonal form in the basis {|00), |10), |01), |11)}, with \ij) = |ion) ® |field). Therefore 

tJit) - cxp|ji^^[: Jo(n):]26ts| 

= e"^°" |00) (00| + e"^"> \ 10) (10| + e">">' |01) (01 1 + e'"^" |11) (11| , (26) 
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where 

000 = 

010 = 

001 = 

011 = nt[Li{rj^)]^, (27) 

being Liijf) the first order Laguerre polynomial and 17 = —g^e^'^'/lS.. In order to have a tt phase gate we need to have 
fit = 2tt and 0oi — 0ii = — tt. This would demand [Li{rj^)]'^ = 1/2, or 77 « 0.54, which is close to realistic values for the 
Lamb-Dicke parameter found in the literature i24ll . 

It is now necessary to give a brief discussion on the experimental implementation of the gate. According to recent experimental 
data the transition from the state Pi/2 to -D5/2 in "'*'Ca+ ions is coupled to the vacuum field with strength g = 2ti x 1.51 
MHz. The losses are mainly due to the cavity at rate k = 27rx41.7KHz and atomic spontaneous emission at rate 7 = 27r x 1.58 
MHz. Choosing A = IO5, the gate operation time is top ~ 9 /^s. We must compare it to the decay rates. Firstly, quite remarkable 
is the fact that the gate works all time with the ion in the ground state. For this reason the proposal is quite insensitive to 
spontaneous emission and this is much desirable. Secondly, when comparing it to k^^ = 3.8 /is, one can see that they are at 
the same order and so that it is still necessary to improve the quality of the actual cavities. However, rather than improving the 
cavity finesse, one could think of using atomic levels which would provide a stronger coupling g. The important point to be 
noticed here is that the robustness against spontaneous emission in our scheme is the key to improve the coupling constant g. We 
remember that this coupling is dependent on the atomic decay rate as g ^ ^ and then one could choose fast decaying excited 
levels leading to stronger couplings and consequently faster gates. It is worthwhile to notice that a phase gate in the Lamb-Dicke 
regime [see Eq. J24H would be much faster than the one we have just discussed. In that case, top ~ 0.62ms which means 
that it is possible to have about 160 cycles per k^I, making our proposals even more interesting for practical implementation. 
Having taken those typical experimental values into account, and considering the robustness against spontaneous emission, we 
find that our proposal opens up new possibilities for feasible quantum computing and also brings the interesting possibility of 
implementing cross-Kerr interactions between bosons in a well developed experimental system. 
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